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Abstract This work is devoted to the exact statistical mechanics treatment of simple inhomogeneous few-body 
systems. The system of two Hard Spheres (HS) confined in a hard spherical pore is systematically analyzed in 
terms of its dimensionality D. The canonical partition function, and the one- and two-body distribution functions 
are analytically evaluated and a scheme of iterative construction of the D -\- 1 system properties is presented. 
We analyse in detail both the effect of high confinement, when particles become caged, and the low density 
limit. Other confinement situations are also studied analytically and several relations between, the two HS in 
a spherical pore, two sticked HS in a spherical pore, and two HS on a spherical surface partition functions are 
traced. These relations make meaningful the limiting caging and low density behavior. Turning to the system 
of two HS in a spherical pore, we also analytically evaluate the pressure tensor. The thermodynamic properties 
of the system are discussed. To accomplish this statement we purposely focus in the overall characteristics of 
the inhomogeneous fluid system, instead of concentrate in the peculiarities of a few body system. Hence, we 
analyse the equation of state, the pressure at the wall, and the fluid-substrate surface tension. The consequences 
of new results about the spherically confined system of two HS in D dimension on the confined many HS system 
are investigated. New constant coefficients involved in the low density limit properties of the open and closed 
system of many HS in a spherical pore are obtained for arbitrary D. The complementary system of many HS 
which surrounds a hard sphere (a cavity inside of a bulk HS system) is also discussed. 



1 Introduction 

The Hard Spheres (HS) and Hard Disks (HD) systems have attracted the interest of many people, because they 
constitute prototypical simple fluids and even solids [3l [5lfT71l29lHilH8l[53] . The extension to arbitrary dimen- 
sionality of this hard spherical particle system has been the object of several studies too [7 lfT4|ll5[ [27 j [3T |[32|l54j . 
Though the apparent simplicity of these systems, only a few exact analytical results for the homogeneous system 
are known concerning mainly the dimensional dependence of the first four virial coefficients |14 [[3mi32j existing 
numerical extensions to higher order [T5ll25] . Exact virial series studies has been also done in inhomogeneous 
systems [8".34'J45j in three dimensions. Recently some attention was dedicated to very small and inhomogeneous 
systems of HS and HD confined in small vessels. The study of small systems constrained to differently shaped 
cavities has enlightening aspects of loss of ergodicity, glass transitions, thermodynamic second law, and some 
other fundamental questions of statistical mechanics and thermodynamic f4j [22|l24|l39|l46[ l4"7 p 55) . 

The exact evaluation of the properties for continuous (non-lattice) inhomogeneous systems of few particles 
is a new trend in statistical mechanics. The analyzed systems are usually HS and HD where the hard potential 
represents the simplest non null interaction and different ensembles approach may be done. An indirect result 
of these calculations is the exact volume, size and even number of particle dependence of low order cluster 
integrals [8l [2n[26lf34U49j . Until now, the systems of two and three equal HDs in a rectangular box has been 
solved [13IJ38! more recently two HS in a box [5lj and two HS and HD in a hollow or spherical cavity [IS] were 
studied. A systematic dimensional approach of the two HS in a box system was also done in j5l], where an 
iterative construction framework was adopted for the increasing dimensional system, but only the two and three 
dimensional systems were explicitly solved. In this work we focus on the analytical evaluation of the Canonical 
Partition Function (CPF) for two HS particles confined in a Hard Wall spherical pore (HWSP) in arbitrary 
dimension D, from now on 2-HS-HWSP. Then this work may be seen as the complement and the dimensional 
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generalization of [l^ . In the rest of the manuscript we will use HS as the dimensional generalization of Hard 
Spheres. With the purpose to avoid any confusion we should mention that pore and cavity are synonyms along 
present work (PW). Sometimes, an empty spherical space inside of a bulk fluid was also called a cavity in the 
literature. 

Prom a more general point of view, the object of this manuscript is the exact analytic study of an inhomoge- 
neous fluid system with spherical interface. This general problem is currently studied because of an incomplete 
understanding of the surface tension behaviour in presence of curved interfaces [9. 11. 12.3 3^43] . It seems that 
the spherical symmetry is the simplest one, and then the principal subject of several works on curved interfaces 
but deviations from sphericity are also studied [33]. The suspended drop on its vapor [OlllOlllQ) . the bubble of 
vapor on its Hquid, the fluid in contact with a spherical convex substrates (or cavity in the liquid) [ 11.12,19] 
[20] . and the fluid conflned in a spherical vessel or pore [TT | fT8 l l40] are different systems in which the spherical 
inhomogeneity of the fluid is central and currently, the study of these systems are converging to the analy- 
sis of the curvature dependence of physical magnitudes [11]. Particularly relevant for PW are such works on 
HS systems [12] and Hard Wall spherical substrates [TT|I20) . PW seen on this context, shows the dimensional 
dependence of an analytical solvable system on this up to date and relevant problem in statistical mechanics 
and thermodynamics. In PW we study a fluid in contact with a hard wall, therefore, we deal with the surface 
tension of a fluid-substrate interface. We know that the point of view of a two-particle fluid may be somewhat 
conflicting. Few-body systems are not gases, nor liquids and neither solids, but we wish to emphasize that the 
toolbox of statistical mechanics must be applicable also to the two particles inhomogeneous system, despite of 
whether it is fluid or not. Naturally, the ergodic characteristics of the system must be considered. We may note 
that the equivalence between different ensembles of statistical mechanics is here invalid. The use of the canonical 
ensemble enable the analysis of such few-body system. In PW it is assumed that the 2-HS-HWSP system is 
well described by the constant temperature ensemble without angular momentum conservation, i.e. the usual 
canonical ensemble, but different assumptions may be done [52]. The macroscopic open inhomogeneous system 
of many HS, interacting with a Hard Wall and even in a HWSP were studied in a virial series way by Bellemans 
[H] and, Stecki and Sokolowski [l5] for D — 3. The virial series or density and curvature power series expansion 
of statistical mechanic magnitudes are of the highest importance, because are a source of exact results, which 
guide the development of the fleld. Therefore, we will make contact between PW and several cluster integrals 
reported in |8jl45j. The flrst terms of the density and curvature expansion of the grand canonical potential, 
surface tension, and adsorption are easily obtained as a by product of the Conflguration Integral (CI) of 2-HS- 
HWSP. Then, we present the value of some integral coefficients related to thermodynamical properties of open 
inhomogeneous systems and comment on this until now unknown dimensional dependence. 

In Section [2] we show how a spherical pore that contains two HS can be treated as another particle. There 
we bring up the central problem solved in PW, the evaluation of CPF and distribution functions of 2-HS- 
HWSP in D dimensions. We also establish a relation between its CI, the CI of three HS and the third pressure 
virial coefficient for the bidisperse homogeneous system of HSs. The exact analytical evaluation of its statistical 
mechanic properties is done in Section [3] where we do a detailed inspection of the highest and lowest density 
limits and make the link with the first terms of the density power series coefficients of some physical properties 
of the bulk fluid system. Principal characteristics of the one body distribution function are also analyzed. The 
Section [4] is devoted to trace the properties of two systems closely related to the 2-HS-HWSP: two slicked 
spheres in a spherical pore and two spheres moving on a spherical surface. The CI of the three systems are 
strongly related as a consequence of the introduced sticky bond transformation. A discussion of the mechanical 
equilibrium condition for spherical inhomogeneous systems and the analytical evaluation of the pressure tensor 
is reported in Section [H The Equation of State, the (substrate-fluid) surface tension and other thermodynamic 
features of the 2-HS-HWSP are studied in Section O where the density and curvature flrst order terms of several 
magnitudes are determined. Several relations with the bulk HS-HWSP open system and its conjugated system, 
when HS particles are outside of a hard spherical substrate, are also provided. Final remarks are given in Section 

m 



2 Partition function and Diagrams 

We are interested in the properties of the CPF of two D-Spherical particles of diameter cr, which are able to move 
inside of a ZJ-spherical pore of radius R' . The total partition function Qd splits on kinetic and conflguration 
space terms. The kinetic term is given by A"^^ where A = (27r/3?i^/m)^/^ is the thermal de Broglie wavelength 
and /3 = (ksT)'^. Then the central problem solved in this work corresponds to the evaluation of CI, i.e., to 
analytically solve Qd 

Qd = QdA^^ ^ J J ^AeseAB dvAdrB , (1) 

CAB = exp[- P(l){r ab)] = 6>(rAB - cr) , (2) 
e, = exp[-f3^jj{n)] = 0{R ~ n) , (3) 



where (f>{r ab) is the hard core potential between both spherical particles, '0(^1 ) is the external hard spherical 
potential, i = A,B and Qd is independent of temperature. The effective radius is R = R' — a/2 and 0{x) is 
the Heaviside unit step function (0(a;) = 1 if 2; > and zero otherwise). Then a variation of a at fixed R does 
not imply a volume variation. For future reference we introduce the Mayer functions / (or /-bond) 

GAB = 1 + /ab , /ab = -0{a - tab) , 



e» = 1 + /^ , ~0{n - a) , (4) 

functions ca, gb and /ab may be seen as overlapping functions, they take the unit value (positive or negative) if 
certain pair of spheres overlaps and become null if the pair of spheres does not overlap, just the opposite apply 
to the /a, /b and cab non overlapping functions. We are interested in the graph representation of Eq. |[T]) then 
we will draw the positive overlap functions, i.e. {ca, cb, — /ab}, as continuous lines, and positive non overlap 
functions, i.e. {—/a, — /b, gab} as dashed lines. The graph of CI is then 



W = 




(5) 



where it is implicitly assumed the integration over the coordinates of both particles A and B. The Eq. ^ shows 
that the system is equivalent to a three particle system. We can apply the in-out relation for three bodies t49j 
performing a simple decomposition over the particle P in Eq. ([5]) 




(6) 



where minus one factors were omitted. All the integrals drawn as not fully connected graphs can be evaluated 
directly because they are factorisable [21]. Focusing on the fully connected graphs, first row relates the config- 
uration integral with the -shape and volume dependent- second cluster integral bniV) of the inhomogeneous 
system, which is also part of the third cluster integral of the homogeneous multidisperse HS system [21]. Second 
row separates the trivial volumetric term from the non trivial area-scahng integral. For only two bodies we have 



= Vd{R) = R''Sd/D, (7) 




-Voia) = -2bD , (8) 

Sd = 2^^/V^(^/2) , (9) 

where Vd{R) and Ad{R) = R^^^Sd are the volume and surface area of the ZJ-sphere of radius R, with Sd the 
soHd angle. The Eq. ([7]) is the accessible volume for a particle in a pore with effective radius R, i.e. the CI for 
the one particle system, and Eq. (JH]) is twice the second cluster integral or second pressure virial coefficient in 
the infinitely homogeneous system of HS which we name bo [21]. With this prescription 5d is a positive defined 
constant. An interesting point is that an inner sphere with radius R — a/2 exist, when this sphere had the same 
radius of the particles, i.e. a/2 or R = a, last integral in first row of Eq. ^ is (minus) 2/3|*'" where f32^^ is 
the second irreducible cluster integral and — 2/3|*''/3 is the third virial series coefficient of the pressure for the 
homogeneous HS system [21] in D dimensions. 



3 The density distribution and CI integration 

With the purpose of evaluate analytically Qd'^g introduce a few geometrical functions. The function Woir, Ri), 
for < r < 2i?i , is the partial overlap volume between two spheres of radius i?i > which centers are separated 
by a distance r, while for — 2i?i < r < 0, it measures the joined volume of partial overlapping spheres. The 
volumes of intersection and union of both spheres are related by WD{r,Ri) = 2VDiRi) — WD{—r,Ri). The 
function Zuir, Ri), is a generalization of the above idea for —00 < r < +00, 



r , if r > 2i?i , 

Znir, i?i) = <^ Wnir, Ri) , if |r| < 2Ri , 
[2VDiRi), ifr<-2i?i. 



(10) 



Similarly, for two spheres with different radii Ri and R2 (assuming Ri > R2) and restricting now to r > we 
have the overlap volume Z]j{r, Ri, R2), 

r , if r > i?i + i?2 , 

ZD{r,Ri,R2)= { WD{r,RuR2),itRi~R2<r<Ri + R2, (11) 
[Vd{R2), ifO<r <i?i-i?2, 

where Woir, Ri, R2) is the volume in the partial overlap configuration. As a consequence of the lens shape of 
the intersecting volumes of two equal and unequal spheres, they are related by 

Woir, RuR2) = \wD{r', Ri) + \wD{r",R2) , (12) 

with r' = r + {R\ — i?2)/r, r" — r ~ {R\ — i?!)/?" and negative values of r' and r" are possible due to Eq. (fTO| . 
We may transform Eq. (fTOl) to a dimensionless function of a; = r/(2i?i) 

r , if x > 1 , 

Zoir.Ri) = VD{Ri)CDix) = VoiRi) X I woix) , if - 1 < a; < 1 , (13) 

[2, if2;<-l, 

Function wd{x) ior < x < 1, measures the overlap or intersection volume between two spheres with unit radii 
separated by a distance x, whereas wd{—x) measures the join volume. Following the same idea Eq. (flTI) may 
be expressed as a function oiy = r/ (2R) (with R — (i?i + i?2)/2) 

r , if y > 1 , 

ZD{r,Ri,R'2) = VD{R)CD{v,A)^VD{R)-^ {wD{y.A), iiA<y<l, (14) 

{wD{A,A),ii 0<y<A, 

and A = (i?i — R2)/(2R). For A = one gets y = x, Cd(?/, 0) = Cd{x) and WD{y,0) = wd{x). Func- 
tion WD{y,A) measures the normalized overlap volume of two spheres with unit mean radius, asymmetry 
A, and centers separated by a distance y with A < y < 1. The analysis of the properties of functions 
{CD{x)XD{y, A),'WD{x),'WD{y, A)} will be left to next subsection 13.11 Now we shall point out that, follow- 
ing a procedure depicted in |49j, we may write down two relevant distribution functions. A pair distribution 
function g(rAB) in which the position of the pore center was integrated, and the one body distribution p{rA) 

m 

ff(r-AB) = Qo'VDiR) eABWDirAB/2R) , (15) 



p{rA) = 2Q-'VDiR) CA 



R + a ( TA R-cT 

( r.n ) 



2R ' ^ \R + (7 R 



(16) 



Function g{r ab) is the probability density distribution of finding both particles separated by a vector tab- 
The normalization equations for these functions are / p{rA)drA = 2 and / g{rAB)drAB — 1- We also introduce 
9{tab) — ff(fAB) Vi)(i?)/2 which is essentially the usual prescription of the pair distribution function in the 
canonical ensemble (see Eq. (29.35) in [21j). Performing the complete integration (of Eq. ifTSj) for example), it 
is found the partition function 



QD=A-^''SDVDiR) / r''-'wDir/2R)dr, (17) 
= A-'^°Vl[R)2^ D r x^-^WD{x)dx = A-^^V§{R)qD{z), (18) 

J z 

where z ~ a/{2R) is positive, qoiz) is the reduced CI, and quiz) = if z > 1. 
3.1 The auxiliary functions wu{x) andwD{y,A) 

The overlap volume between two spheres with unit radii in any dimension, wd{x) for < a; < 1, introduced in 
Eq. m is [7] 

WDix)^h-^,2{{D + l)/2,l/2), (19) 

being Ix{a, h) = Bx{a, b)/B{a, b) the normaHzed incomplete beta function, defined in terms of the beta function 
B{a,b) = r{a)r{b)/r{a + b) and the incomplete beta function B^{a,b) = dtf-^il - tf-^ [2\. For future 
reference we introduce the shortcut Bd = l/i?((il'-f l)/2, 1/2) = iS'd+i/S'd+2. We extend definition lfT9|) following 
Eqs. (fTOlfTS]) to w_D(a;) — 1 = —{wd{—x) — 1) for — 1 < a; < 0. In turn, due to the relation /i_j.2(a, b) = 1—Ix2{b, a), 
for — 1 < X < 1 we have 



wd{x) - 1 = - sign{x) 42(1/2, {D + l)/2) , 



(20) 



where sign[x) = 1 if x > and sign[x) = — 1 if a; < 0. Therefore, for £> > 0, woix) is analytic for — 1 < a; < 1 
and both {wd{x) — 1, C,d{x) — 1} are odd functions. It is also possible to construct wd{x) from the recurrence 
relations f7| 



w-i{x) = 1 , 
2 

wo{x) = — arccos(a;) , 



wd{x) = wd-2{x) - x{1~ x^) 



2\iD-l)/2 



B_i = 0, Bo = I/tt and Bjj = D/ (2ttB]j^i). Expressions for D ~ 1, 2, 3, 4, and 5 are 

wi{x) = 1 — X , 

^2(2;) — — arccos(x) — a; (l — a;^)^^^ 



W4{x) 



/ ^ 1 3 1 
W3[X) = ^^2^+2' 

5 2 



arccos(a;) 



W5{x) 



15 



1 - —X + -x-' 



t:X 



(21) 
(22) 

(23) 

(24) 
(25) 
(26) 
(27) 
(28) 



8 " ' 4^' 

As far as we will need the asymptotic analysis of the wd{x) function, it is resumed here. In the a; ^ and 
X — > 1 limits we have respectively 



wd{x) = 1 - 2Bz3.e + ^ ^ Bdx^ - 0^{x) , 



wd{x) 



2(^+3)/2Bi3 

D + 1 



(l-x)(^+i)/2 (l + 0(l-a;)) . 



Following Eqs. lfT0|) - (fT3|l woiy,^) may be written as 



(29) 
(30) 

(31) 



with x' = r'/2Ri = cos(0i) = (y^ + A)/[y (1 + A)] and x" = r"/2i?2 = cos(6'2) = (y^ - A)/[y (1 - A)], here Oi 
is the angle opposite to the side of length Ri in the triangle {Ri,R2,r), i.e. cos(0i) = (— i?f + i?| + r^)/(2i?2r) 
and cos(6'2) — (-Rf — i?2 + ^'^)/(2i?i?'). A recurrence relation for wd(j/,^) is derived in the Appendix [Aj From 
Eqs. (|2T1) - I|28p or from Eqs. (|132p - i|139p of the Appendix [A] we have that the first functions of the series are 



w^,iy,A) = 1/(1 -A^), 



wo(y, ^) = vo{y,A) 
1 



voiy, A) = -{I + Z\)^ arccos( ^ 

TT y(l + A) 



y 



(l-Zi)^arccos(- 
71" 2/(1 --4) 



leading to the following expressions for D = 1, 2, 3, 4, and 5 

wi{y,A) = 1 - y , 



w2{y,A) 



{{l-y^){y^-A^)Y" + V2(y,A), 



W4{y,A) 



W3{y, A) = (2y)-i(l - y)^2y + y^ - 3A^) , 
-^y-' ((1 - y'){y' - A^)f' (2y^ - 5y'{l + A') 4A' 



vi{y, A) 



W5{y, A) = (2y)-^(l - yfiSy^ + 9y^ + 3y^ - 10AY{3 + y) + 5A^il + 3y)) , 
where Eqs. (I32l)-ll39l) apply for max{0, A) <y <l. 



(32) 
(33) 
(34) 

(35) 

(36) 
(37) 
(38) 
(39) 



Figure 1 (color online) Reduced partition function qoiz) as a function of 2 = a/2R for several dimensions. From left 
to right D = 1, 2, 3, 4, 5, 10, and 50 last two in dashed line. 



3.2 The reduced CI 

As shown in the Appendix [HI the reduced CI qoiz) defined in Eq. (fTTI may be written as 

qniz) = ud{z) - {2z)^wd{z) , 
for < z < 1 while qoiz > 1) = (7d(1) = 0, with 

uniz) = /i_,2 {{D + l)/2, {D + l)/2) , 
valid for D > but not for D = —1. Quantity uu may be derived from the recurrence relation 

= 1/2, 

uo{z) ~ (2/7r) arccos(z) , 

uoiz) - UD-2{z) + z^'-'il - z^i^'-^^l^i}- - ■?)r{D)/r^{{D + l)/2) . 



The CI for LI = -1 and are 



<l^.iz)^\il-\) 



qo{z) = , 



(40) 

(41) 

(42) 
(43) 
(44) 

(45) 
(46) 



while for D = 1, 2, 3, 4, and 5 are 



qi{z) = (1 - zf , 



(47) 



2 r 

92(2) = - (1 - (2zf) arccos(z) + z ^1 - z^ (l + 2z^) 
q3{z)^l-{2zf + 9z^-2z\ 



94 (^) 



(1 - (2z)^) arccos(z) + -z yj\ - z2 (3 + 2z^ + SGz"* - 16z^) 
3 



g5(^) = 1 - {2zf + 50zS - 25z* + 6z 
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(48) 
(49) 
(50) 
(51) 



It is apparent that CI for odd D is polynomial with order 2D and integer non null coefficients at terms of order 
0, D and 13 + 1 + 2fc with /c = 0, 1, 2, (D — l)/2. However, for even D partition function is not polynomial. 
These and other interesting properties may be derived from the series representation of the incomplete beta 
function (from [l] and Eq. (26.5.4) of [2j). Function qj:i{z) is plotted in Fig. [T] for several values of D. For all 
D > we have 50(0) = 1, 9r>(l) = 0, and for large D limit qo^ao^z) = 0(2"^/^ — z). Let us now look at some 
asymptotic behavior. 



Large cavity limit: The infinitely dilution or large cavity limit z ^ (i.e. a/R ^ 0) can be studied from the 
series representation of the reduced partition function qoiz) near z ~ 0, 

oo 

qoiz) = 1 - z^C,,o + z^+'Cs + z°^^C,^2 + E ^''+'+''C'i,2+fc , (52) 

fc=i 

where we introduce a set of dimensional dependent constants {C} that will denote some power series coefficients 
through PW. Here, Ci^o = 2^, Ci,i =2^Bd^ and Ci,2 = -2^Bd^^§j:^. For odd £> Eq. JH) is effectively 
an order 2D polynomial. As we will see later, the constant coefficients {Ci^o, C*!,!' ^1,2} are involved in physical 
properties of the equivalent many body system. From Eqs. H]) and (fTSl we obtain 

Qd = V^{R) - Vd{R) 2bD + Ad{R) 2aD - JJd{R) 24'^ + Od-^{R) , (53) 

written in terms of the extensive squared mean curvature JJd[R) = Ad{R) ■ [D — 1)^ • R^^ . The coefficients of 
area and curvature are 

=foD+i(2^)-\ (54) 

4'^-fei.+3§ij(2^3.^)-\ (55) 

and (J^"* — corresponds to the absent term proportional to A£){R) R^^ . We may remark that curvature 
dependence is neither proportional to total nor to Gaussian curvatures, j = {D — 1)/R and k = i?^^^^^) 
respectively. As may be expected VD{R)2bD is the first non ideal correction, and first sign of inhomogeneity 
and curvature dependence appears in the next two terms. They are deeply connected with 6_D+j the second 
cluster integral in higher dimensionality. Therefore, we are showing a direct relation between the -intrinsically 
inhomogeneous- properties of 2-HS-HWSP and the low density limit properties of HS homogeneous system 
in higher dimension. Equations l(53l [5i|) and l(55|l also mixes the properties of systems with odd and even 
dimensions. Finally, note that Eq. ((53| is exact for Z? = 3 without the order D ~ 5 term [49j and coincidentally 
~ k with the constant value JJ3{R) = 2^'k. Further, JJd{.R) ^3xId=3 = ■U?,{R) + '2K-i{R) with the extensive 
gaussian curvature Kd{R) — A]j{R) ■ k. Fixing R = a the inner sphere and both particles have the same size 
and we obtain from Eqs. ^ [53|) 

2/3r \d=3= -(263)' + Asia) 2a3 - J73(a) 2(5^ , (56) 

where the second irreducible cluster integral /Jj*'^ which involves three nodes [2T|, is now written in terms of 
geometrical measurements of the cavity's boundary {As{a), JJ^ia)} in three dimensions and two body integrals 
&£) in D > 3. Baus and Colot [7] found (32^^ \d= — (&_D)^t«D(l/2). These type of relations may be interesting for 
higher order integrals. The three body integral (two HS plus the cavity) Qd and its moments are linked with 
several three body integrals describing physical properties of the HS inhomogeneous fiuid in the low density 
limit, inside a spherical cavity and even outside a spherical substrate, which will be discussed later in Sec.[6l 



Small cavity limit: Looking at the opposite situation of high density or caging limit we found the final solid, 
i.e., the densest available configuration of the system. This limiting behaviour has been extensively studied for 
the HS homogeneous solid system, but also, in small systems |39P49). The caging limit of 2-HS-HWSP is obtained 
at the root z = 1. The reduced partition function has an interesting series representation in the neighboring of 
z = 1 (valid for z < 1) 

quiz) = (1 - z2)P+3)/2c2,o (1 + C2,i (1 - z^) + 02(1 - z^)) , (57) 

with 6*2,0 = Ci_i/{D + 3), 6*2,1 = 2(0+5) ■ shows that CI goes to zero as (1 — 2)(^+'^)/2 when system becomes 
caged. The partition function is then 

Qd ^ a3(^-i)/2(i? - (7/2)(^+3)/2^3,o(l - f7"'C3,i(i? - cj/2) + 62(1 - 2i?/a)) , (58) 

here 6*3,0 = C2fi2^-^ {Sd / Df and 6*3,1 = 2-i(-9 + 8 62,1 + 5£i). Since Qd is identically zero for R < a/2 Eq. 
(|58|) impHes that is non-analytic at R~ a/2. For odd D, the derivative of order 2 -I- (D -I- l)/2 and beyond 
are zero. However, for even D the derivatives of order 2 + D/2 and bigger involves an infinite discontinuity. 
When a closed system of hard particles is caged as a consequence of the high confinement, the spatial degrees of 
freedom that becomes lost are related to a zero measure set in the CI integral similar to ([T]). The consequence 
is that Q becomes zero, and the signature of the frozen spatial freedom is the order of this root. We introduce 
the number of spatial lost degrees of freedom LDF and the complementary number of kept degrees of freedom 
KDF, their addition provides the total spatial degrees of freedom LDF + KDF — N ■ D where iV = 2 is the 



R/0 

Figure 2 (color online) Characteristic values of density profile. Densities, at the central plateau po (continuous line) and 
at contact pc (dashed line) as a function of pore size. From top to bottom dimensions D = 3, 5, and 7. Rough density p 
for D = 3 is showed in dotted line. 

actual number of particles. We relate LDF with the exponent {D + 3)/2 in ((58)) . For the studied system of 
2-HS-HWSP in the caging limit we obtain 

LDF ={D + i)/2 = l + {D + l)/2, (59) 

KDF = 2D - LDF ={D- 1)3/2. (60) 
In next sections we will study LDF and KDF in other situations. 



3.3 The one body distribution 



The one body distribution function are qualitatively similar for all the dimensions and was previously described 
in detail for D = 2 and 3 [49] (Fig. 2 therein). A peculiarity of p{r) is the plateau of constant density that 
appears for R > a, and which extends from the center to r = R — a. For a/2 < R < a a null density plateau 
develops in the range < r < a — R. Therefore, the central plateau density is 

Po = 2Q^'Vd [1 - {2zf] ^2{Vd~ 2bD) /Qd , (61) 

if i? > cr, while po — ii a/2 < R < a. Other interesting magnitude is the density at the wall or contact density 

Pc = p{r = R), 

Pc^2VdQd' l^^-^[wD{l-2z^) + {2zfwD{z)]^ ■ (62) 

These quantities are shown for several dimensions as a function of pore size in Fig. [2l We introduce the rough 
or mean number density p — 2/Vd{R) for comparison. We would like to emphasize that p[r) is a discontinuous 
function at the cavity surface falling down to zero outside the cavity and, also, is non-analytic at r = Abs{R — a). 
The function p(r) a regularized version of p{r) at r — R, may be introduced by the short-hand 

p{r) ^ eA{r) g{r) , (63) 

which must be understood with the help of Eq. lfT6|) . It is similar to p(r) but includes its analytic continuation 
for R < r < R + a and therefore is a smooth function at the hard wall ([16], p. 166 therein, this procedure was 
also used in the case of homogeneous fluids, see Eq. (6.11) in |5j). 



4 Other Closely Related Systems 



It is possible to establish interesting relations between the CPF of 2-HS-HWSP and that of other systems. One 
of them involves the one stick or dumbbell in the spherical pore, where we think the stick as a rigid body formed 
by two sticked HS. Other system is that composed by two hard spherical bodies which are constrained to move 
between the HWSP and an inner hard spherical core. For sufficiently large hard core radius we find the limiting 
case of 2-HS that are able to move on the surface of a sphere. We may imagine that this system is constituted by 
two HS sticked to the surface of the HWSP, all immersed in a D dimensional euclidean space. Therefore, it has 
an effectively reduced dimension of D—1. Both systems may be seen as originated in the 2-HS-HWSP and they 
emerge as a consequence of the addition of a new sticky property. Interestingly, this property is easily introduced 
by a simple transformation. In PW we will not make a general digression about this transformation, but we 
simple state that a sticky-bond transformation is performed when one or several of the e or f functions in ID 



[TSl) is transformed in a Dirac delta function or s-bond. In a system of many hard bodies in a hard pore a sticky 
transformation may be done for example by sticking one particle on the surface of the pore, or sticking particles 
together. The sticky-bond transformation is interesting because in an iV-body system the most compact spatial 
configuration, or final solid may be obtained through a series of such operations. It may also mimic a nucleation 
or condensation process, the adsorption on a surface and even some chemical reactions. In this section the 
temperature and kinetic factor are not relevant, then we simply make A = 1 and then CI and CPF are equal. 
We may mention that the sticky bond implies that the two bodies must fix their distance, then, even when the 
concept is related with the adhesive hard sphere bond of Baxter [6] they are not equivalent. 



4-1 The one stick in a HWSP Partition Function 

When two HS become stickled together they conform a rigid body because the s-bond fix the inter-particle 
distance to a, they form a stick. The sticky-bond transformation may be done by applying the derivative da to 
Qd on Eqs. CD [111) 

Qo/b^CbdaQD, (64) 

the 6-subindex label means body-particle opposite to point-particle which may be used to describe a HS. The 
constant Cb = found from the large size limit, obtaining 

Qo/b^ SdQdM^^) = SdVd{R)wd{z) . (65) 

Iterative construction of wd{z) and limits z — > and z ^ 1 were described in Section ISTTl Limiting behaviour 
of the partition function Qo/t is, from Eq. (|29l) 

D + 1 

= Sd {Vd{R) - Ad(R)<tBd/D - JJd{R) a^Bz, {A\D{D - 1))-^) + Od-5{R) , (67) 



Qo/t - SoVoiR) - AoiR) 5i5+ia(27r)-i - JJd{R) So+sj^^H^^^^^r' + Od-:.{R) , (66) 



Note that Vd{R) is the volume for one HS particle in the spherical cavity but Vol = Qo/b/SD is the actual 
volume of the pore for the rigid particle composed by two point particles with fix separation cr, they are strongly 
different only for small R in fact 



Vd{R) - VoljR) _ a_ 



R) 4! 



Od^). (68) 



The first correction to the one HS CPF is related to the Volume coefficient of Q^i and the next one relates to 
the Area coefficient as a consequence of Eqs. l(53t [65 |) . The relation may be compactly written o-°~^+™S'D+m = 
D-\-m and is a direct consequence of Eqs. ((TlIHl)- The power at which Q^i/b .i.e LDF is now, from Eq. 

ISQD, 

LDF^{D + l)/2. (69) 

When this system becomes caged it conforms a linear rigid rotator. The difference between LDF in Eq. l(69|) and 
LDF in Eq. (f59)) is attributable to the relative distance between both HS. Then we interpret that LDF oi Qd 
corresponds to one degree of freedom from the relative distance between both HS and {D + l)/2 corresponding 
to the LDF of Qo/b related to the rigid rotator degrees of freedom. 



4-2 The spherical pore with a hard core 

Other closely related system is that of two particle confined in a spherical pore with a fixed central hard core. 
The main interest to solve this pore shape is the study of the dimensional crossover between D and D — \ 
dimension, which happens when the internal core becomes so large that the two HS may only move on the 
surface of the D-sphere, an effective D — 1 curved dimensional space. This type of non planar surfaces was 
introduced by Kratky |26j in the study of the HS homogeneous systems. Exact results for two, three and four 
hard particles confined in this non eucHdean surface embedded in £> = 3 have been reported [41]. Besides, 
the 2-HS-HWSP with a hard core (2-HS-HWSP+HC) is by itself interesting cause it is a highly non trivial 
concave pore. The actual volume of the system is Vd{R) — Vd{R — h), being R — h > the radius of the 
core and < h < R. The CI of 2-HS-HWSP+HC was written in Eq. HJ but Eq. © must be replaced by 
Bi = exp[—f3ip{ri)] ~ 0{R — ri) — 0{R — h — r.^). Expanding the product of Heaviside functions we obtain 

p2R. l-2(R-h) 

QD{R,h) ^ SdVd{R) r°-^WD{r/2R)dr + SoVoiR-h) r^-^WD{.r /2{R ^ h)) dr 

Ja Jmin{cr,2{R~h)) 

-2SDVD{R-h/2) / r''-^i:D{r/{2R-h),h/{2R-h))dr, (70) 

J min{a,2B.— h) 



where we have made expHcit the conditions on the integration intervals. From them, we observe that CI breaks 

in three branches depending on cr < 2{R — h) or 2{R — h) < a < 2R — h or 2R — h < a < 2R. Last row, 
also shows that CI separates in two branches according to Eq. (fT4| depending on cr < /i or not. Equation i[70|) 
becomes 

QoiR, h) = Vl{R) qoiz) + V^{R- h) 0{2{R ~ h) ~ a) W2(i? - h)) 

-2Vl{R)e{2R^h-a)2''D J y""-^ woiy, A) dy (71) 

-26)(2i? -h- a)e{h - a) Vd{R - h){VDih) - Voicr)) , 

being R ^ R — h/2 and r = max{A,a/2R). It is clear from Eqs. ((7T| [3T|) that at each branch Qd is a finite 
power series in R, h and a only for odd dimension. The dimensional systematic of the integral in l|7ip becomes 
a complex task, then we only analyse with more detail D = 3 and some characteristics for D = 2,4:, and 5. For 
D = 3 we evaluate the integral to obtain an expHcit expression for the partition function, here we present the 
result for a>2{R~h), 

' Vof ~ Vol 263 + Ar 2a3 ~ Cur 25^^ , ifh>a, 

QsiR, h) = { ^^^2 _ ^ ^^^^ (8 - 6A + A3)/3+ ^^^^ 

-Citr24^^A3(8-9A + 2A3) , iih<a, 

where A = h/a, and Vol, Ar and Cur are the measurements of volume, surface area and boundary quadratic 
curvature for the actual pore, i.e. Vol = ^3(7?) -Vz{R- h), Ar = ^3(i?) + A-i{R- h) and Cur = JJaiR) + 
JJz{R ~ h) — 2^Tr. Noticeable, the partition function is polynomial in each domain, and has continuous second 
derivative in ft. = cr but a discontinuous third derivative. It is surprising that the volume coefficient 2&3 and 
the entire partition function when h > a looks exactly equal to Q3 (see Eq. l(53|) and Ref. (1^) with different 
Volume, Area, and Curvature measures. We also have established that the central position of the fixed hard 
core is not essential, first row of Eq. (|72|) is valid for any internal fixed hard core as long as its boundary is 
separated from the outer spherical pore wall a distance > cr. The limiting behaviour as ft ^ to first non null 
order in ft is, from Eq. (|72l) 



QsiR, ft « a) = ft2 (AliR) - A^iR) na^) + O^ih) = h-^S^AsiR) {R^ - {a/2f) + O^ih) , (73) 

coincidentally 62 ~ Trcr^ is equal to the first non ideal gas correction. Here Q goes to zero as ft^, i.e., the system 
lost two degrees of freedom. It is an expected property, which must be valid independently of the dimension 
and even the number of particles, Qd{R) — ft^ for R large enough. After h goes to zero, Q/h^ goes to zero as 
(i? — cr/2) it is the caging limit. In the 2-HS-HWSP in Z? = 3, three degrees of freedom are lost when R cr/2 
(see Eqs. (|49t [57 t [SS)) ) as the system becomes caged, in the spherical pore with an internal core they are lost in 
two steps. 

The study of the systematic dependence on the dimension number of the ft ^ limit is not easily available 
from Eq. ((70| . However, from Eq. l[7T|) we find that the partition function for h is Qd — h^Qo/s + 03(h) 
with Qd/s the configuration integral for both HS confined to the surface of the £)-dimensional sphere also known 
as the calottes problem [H]. By means of two stick transformation, the Qd/s may be expressed as 



Qd/s = Sd dR,idR2Q d{Ri, R2)\ri^r2=r = 9fli9fl2 / Woir, i?i, i?2)r^~^ dr 

J a 



, (74) 

R1=R2=R 



where Qd{Ri, R2) is the CI in Eq. (fTSj) but with a different confinement radius for each particle, which is known 
for £> — 2,3 [49]. Here, we transform the ca, cb Heaviside functions on Eqs. |(ll[l7l) into Dirac delta functions 
through the derivatives —dm, sticking both HS on the surface. For practical purposes is much easier to evaluate 
Qd/s from 



Qd/s = Ad{R) Sd-iR""-^ I sin^-2(^) ^ ^^g^ 

Jo 

QD/s^AliR)uD-2{z), (76) 

where Ad{R) is the system volume and the minimum angular distance between particles is d = 2arcsin(z). 
Here ud-2{z) is the reduced partition function which plays the same role that qd played in the HWSP, then we 
define qD/s{z) = ud-2{z). As far as the iterative construction of the ud{z) functions was previously described 
we present the function for 13 = 1, 2, 3, 4, and 5 



qi/s{z) = 1 , 



(77) 



(l2/s{z) = — arccos(z) 

TT 

q?,/s{z) = 1- , 



qi/s{z) = 



arccos' 



(z) + Vl - z'^z{l - 2z^) 



(78) 
(79) 
(80) 
(81) 



an extra 1/2 factor must be considered in Eq. ([77]) due to a kind of ergodic to non-ergodic transition. The 
asymptotic behavior may be accomphshed with the approximate series representations for qjj [l] 



qD/s{z) ^1- z 



D-l 



2(D + 1) 



23(1? + 3) 



(82) 



qn,s{z)^{l-z^f^-^y^C, 

where C4 = ^i^^T^^ "^^2"'")] ^ ^^"-^ truncated terms in Eq. ((82|) are even powers in z. A noticeable charac- 

teristic of both series is that they have the same numerical coefficients and both are polynomial with 2{D — 2) 
degree for odd dimension. The approximate series for Q/j are then 

Qd/s = Al{R) - Ad{R) 2bD-i + JJd{R) 26d+i §3^(2^'^)"' - Od-5{R) , (84) 

where we note that bo-i — 2'KaD-2 and the curvature coefRcient is 67r(5j^^2 in concordance with Eqs. (|54l ISS)) . 
and 

where C5 = 2^(^^^'B£)/I? and the low density limit implies a vanishing curvature limit too. From Eq. fM)) we 
may obtain the first cluster integral correction due to the space curvature for this non euclidean container or 
spherical boundary conditions. We may compare the first cluster integral in the {D — l)-euclidean space 2bD-i 
with the cluster integral in the surface of a sphere in 15-dimensions 2hQ/s = {Qd/s ~ A\){R)) / ■^Di.R) ^ only 
for D = 3 we have 2b£)/s — 26d_i, for any dimension we obtain lim 2b£)/s = 26d_i. Partition function for 

D = 2, 4, and 5 are 



(i^-3)(i^-l) 
2(15 + 1) 



(I - z') + 02(1 ^ z^) 



(83) 



1 - CT 



_, {D~1){7D~9) 
2{D + 1) 



{R- a/2) + a-'^OiR- a/2f 



(85) 



2R 

TT 



h/s = Al (R) - A2 (R) S2 ^ arcsin(z) , (86) 



)4/5 = (7ri?3)2(i „ arcsin(0) + z(l - 2z^){l - z'^fl'^) , (87) 



Qs/. = AjiR) - A.,iR) 71^4/2 + A (i?).i?- W5/12 

= ^^V2y R^iR-^fiR+^fi2R' + a'). (88) 

Analyzing the way in which Qd/s goes to zero we obtain from Eq. (|83l) that, 

LDF= (D-l)/2= (D + 3)/2-2, (89) 

i.e. two LDF from Qd corresponds to the confinement of both particles on the cavity surface. Now, we are able 
to extend the results about LDF in Qd, Qd/i and Qd/s we may argue that if we stick both particles between 
them and to the surface we find 

LDF= (D-3)/2== (D + 3)/2-3, (90) 

where each stick-bond transformation has reduced the LDF on one unit. The idea is that the same final solid- like 
configuration can be obtained in several ways, but its state properties can not depend on the particular taken 
path. Therefore, to stick particles between them, next stick each particle to the surface and finally cage the 
system must produce the same CI (basically the system free energy) that is obtained directly from the caging 
of the 2-HS-HWSP. 
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Table 1 Several constant coefficients of the series representation of CIs and other functions. Some values for D = — 1 
and correspond to a suitable limiting behavior. 



5 Mechanical equilibrium and Pressure tensor 

To make a more complete and microscopic characterization of the system it is necessary to evaluate the pressure 
tensor. Briefly, pressure tensor, density and external potential are related as a consequence of the mechanical 
equilibrium by 

V • P + pV^ = , (91) 

where the exphcit position dependence of the magnitudes has been dropped. In any inhomogeneous system the 
tensor can be split into 

P = /?-Vld + P^, (92) 

where Id is the DxD identity matrix and P*^ is the interaction part of the tensor [36j- In systems with spherical 
symmetry P'^ is diagonal with only two different components [9}l40] the normal and tangential components 
and 

P^ = P^rr + P^{v^ + hOi + ... + 00-200-2) , (93) 



where (p, Oi, ..., 0o-2 are the angular versors. Then, for systems with spherical symmetry Eq. I|9ip may be 
written as 

drPN + ^^—^ {Pn - Pt) + pdri^ = . (94) 
r 

Further simpliflcations apply to a hard wall container, there p(r) is discontinuous at r = i? but g{r) is not (see 
Eq. we obtain from Eq. (|94)l 

drPN + ^^{PN-PT) = -r'S{r-R)g, (95) 
r 

drP^ + ^^—l{P^~P^) = -p-^e^{r)drQ, (96) 



both equations are equivalent to Eq. l(94|) and then necessary conditions for an acceptable pressure tensor 
deflnition. Nevertheless they are formally strongly different. The inhomogeneous term in the differential equation 
((95|) shows at the boundary a divergent singular contribution to the total pressure components with zero 
contribution from positions inside and outside the boundary. Although, inhomogeneous term in Eq. I|96p shows 
at the boundary a discontinuous singular contribution to P^ components, with non zero contribution from 
inside points. We will regress to this issue in a forthcoming paragraph. 

Turning now to the interaction part of pressure P^, is known that its detailed expression is non unique. 
Particularly, different possible definitions of P^ produce different values of pressure tensor in inhomogeneous 
fiuids [9]. In this work we adopt a pressure tensor definition extensively utilized in MD simulations |36], the 
components of the pressure tensor for the two body system are 



Pi;,{r)^{rl^Fi^{v^^)5{v-vp^)) , 



(97) 



been r" = r • a and F^g (r) — — ^ -77 the component of the force between particles in the 6 direction (see Eq. 
(3.20) in Ref. [36]). It has been argued that Eq. (f97| can be obtained from the Irving Kirkwood pressure tensor 
with the assumption of short range interaction [36] . Interestingly, the adopted definition for the pressure is an 
intrinsic two bodies emergent property which depends only on the position of two particle coordinates and does 
not depend on any choice of an integration path. In the present problem is important to notice that position 
of both bodies are the point of pressure evaluation r = ta and the integrated position tb in the mean value 
of Eq. ((97|) . It is a desirable property of P'-^(r) in any few body system that the microscopic configurations 
with zero probability to find a particle in position r must not contribute to the pressure. This property is not 
accomplished by the Irving-Kirkwood choice for the pressure tensor. Equation ([QT)) may be written as 



Pl{r) = (/3Qz5)^^eA(r) 



eABee 



AB' AB 
kABl 



(98) 



where tab = r — tb and Bab ~ <5(rAB — o-). Changing the integration variable to tab and integrating over the 
radial coordinate one gets 



BA (r) J esirB) ^ab ^ab c^^ab , 



with tb = |r — crfABi and = f • b. The two independent components of this tensor are 

^wW = (/3Qi?)"'eA(r)a^5D-i J cb cos'' 9 sin^-H d9 , 

P¥ir) = {l3QD)-'eA{r)a'' Sd-i{D - ly^ /eBsin^^rf^, 



(99) 



(100) 



(101) 



cr, 



where the integration interval (0, tt) may be reduced to (6'mi„, tt) by the effect of cb, with 9min = if r < i? 
9min = TT if \r — a\ > R, and 9min — TT — arccos(— i?^ + ct^ + r^)/(2rcr) otherwise. Both integrals expressed in 
terms of known functions are 



P}^ir) = {PQDy'eA{r)2-WDia) [Z? WD-2(cosa) - (I? - 1) ^^(cosa)] , 



(102) 

P¥{r) - iPQur^eAir) 2-Vi3(a) woicosa) , (103) 

where cos a — 1 if r > R + a, cos a = —sign{R — cr) if r < R — a and cos a = (— i?^ + cr' + r')/{2ra) if 
R — a < r < R + cr. Figure [3] shows both pressure tensor components as a function of position into the cavity. 
We observe that when the pore is big enough, R > a, a pressure plateau develops at the center of the pore in 
the range < r < i? — cr. In the region R — a < r < R the inhomogeneous pressure region develops. The shape 
of pressure tensor and density profiles are simply correlated (see Fig. 2 in [49]), their constant value plateau 
and inhomogeneous regions coincide, even more, when plateau density becomes null, pressure goes to zero too. 
In the homogeneous region the normal and tangential components of the pressure tensor become equal, and the 




r/c 



Figure 3 (color online) Pressure tensor components. Continuous line corresponds to Pjv and dashed line to Pt- Two 
cavity sizes R/a — 0.9 and R/a = 1.2, and two different dimension D are shown. The arrows indicate R/a values. 



pressure tensor reduces to a constant scalar pressure Pq [9]. In this case Eqs. I|1021 IIOS)) lead to 



(/3Qd)"'26b 



finally, according to Eq. ([921) the pressure in the homogeneous density plateau is 

PPo^PO + ^^ 2(Vd - bD)/QD , 



(104) 



(105) 



for R> a while (3Pq — for a /2 < R < u . The value of the pressure tensor at contact can be obtained from Eqs. 
(|102[ llOSl with r = R and cos a = z = a/{2R) but no further simplification may be done. In Fig.|4]we plotted 
the characteristic values of the pressure tensor. There we can observe that the maximum of Pq is attained at 
smaller R than the maximum of the density (see Fig. ^ . We have verified that in the inhomogeneous region 
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Figure 4 (color online) Characteristic values of pressure tensor. Pressure, at the central plateau Po (continuous line) 
and contact value (r = R) of both components Pat and Pt (dashed line and dash-dotted line) as a function of pore size. 
From top to bottom dimensions D = 3, 5, and 7 in triangles, squares and circles. 



abs{R~a) < r < R, Eq. l(96l) with Pn,t from Eqs. I|102lll03l) is false. Two probable reasons may be argued, the 
invalid short range hypothesis for the HS potential for nonuniform density, and / or the incorrect pressure tensor 
definition of Eq. (|97l) due to the hard wall boundary condition (see Eq.(3.5) in |,40j). As far as the mechanical 
equilibrium of Eq. l(96|) is still valid in the homogeneous plateau, and pressure should not strongly depend on the 
definition details in this region, we accept the validity of the obtained pressure tensor for < r < abs{R — a). It 
is particularly interesting to analyse the point r = R where we find drQ > 0. We may note that the discontinuous 
behaviour of Pj^ (but not of Pj^) at r = R violates Eq. ((96l) . The existence of the Heaviside factor e^(r) on 
the right hand side of Eq. l(96|) must be balanced with the same global factor at the left side, and forbids the 
appearance of an uncompensated Dirac delta. However, the derivative of a discontinuous P^ just produce a 
singular Dirac delta at r = i?. If we assume the validity of this equation one obtains that the discontinuity is 
completely unphysical. Thus, we need the continuity of the interaction normal component of the pressure tensor 
to overcome the mismatch. 



6 The Equation of state and the Laplace Equation 

We are interested in the equation of state (EOS) of the 2-HS-HWSP, i.e. the thermodynamic description of the 
complete inhomogeneous system. We adopt a point of view usually taken in spherical droplets which makes a 
description in terms of the properties in the homogeneous regions of the system. In this section we drop any 
explicit unnecessary D subindex and dependence on cavity size R, then V = Vd{R), Q = Qd{R) and so on. 
The properties in the homogeneous region may be found from Eqs. l(6T | llOSp . From these Eqs., for R > a we 
obtain 

BPq b p , , 

Po 2 I - bp 

= l + \{V/2b-l)-\ (107) 

where the compressibility factor at the constant density plateau is a simple function of p or V. We may note that 
both expressions diverge aX R = a and V — 2b. For this pore size the plateau of constant density vanishes and 
Po = 0. Interestingly, in Eq. (|107p we have found the bulk van der Waals EOS without the attractive term, the 
only fiavor remaining the two body system is a 1/2 factor. Owing to Eq. I|106p does not depend on geometrical 
parameters which resembles the cavity's spherical symmetry it should be valid also for 2-HS confined in pores 
with other geometries. In Eqs. I|1061 1107]) Pq, pq and p (or V) are not independent variables, even, they are 
functions of system's size R. The Eq. (|106p can be written in two forms which resemble the contact theorem for 
the bulk homogeneous HS system (see Eq. (2.5.26) in [TBj and also [42] ) 

pPo = po + bp^ 4^ , (108) 

W[Z) 

b o q(z) 

= ^"+2^0(1%^' (10^^ 




Figure 5 (color online) Pressure in the homogeneous plateau. Prom top to bottom D = 2,3, 4, and 5. Two different 
choices of density parameter are depicted. At left hand side (a), in dashed line we adopt the po as density parameter. At 
right hand side (b), in continuous line we adopt the rough density p. 



where in fact g{a)/w{z) — {2q{z))~^ and then g{(T)/w{z) — > 1/2 when oo. The Eqs. I|1081I109|) make sense 
only for < z < 1/2 or cr < i? < cx), where for the smallest pore {R = a and z = 1/2) po — 0. In such a case, the 
term which multiply in Eq. I|108p takes a finite value, but, the term which multiply pg in Eq. I|109p diverges. 
However, Eq. (|109p goes to zero. This shows an only apparent different behaviour, because Eqs. i|1081 11091 are 
different representations of the same equation. At large R the overall and plateau densities are related by 

Po^p-^p'A{R)^p~CeP^+'^'' , (110) 

with Ce = {D^-^2-^SDY'"aD and A{R) - p^^'^^^^''^ - p-(D-i)/D j^^. £) > Polynomial expressions for the 
pressure as a function of density, may be found by considering the firsts terms of Po as a density power series 

PP^^ p^+^-pl+^^f^,+ ^-pt{2h' + aA{R)) , (111) 
:^P+\p^ + \p'{b^-aA{R)) . (112) 

Here we choose two density parametrization, the plateau density and rough density, po and p respectively. The 
expansions in Eqs. I|1111I112|) end at the order of the first signature of the inhomogeneity. Both equations show 
that the first correction to the ideal gas behaviour is positive and involves the expected closed system correction 
[23]. Besides, the Eqs. (|llll 1112]) involve terms with higher powers in the density than two. This feature may 
sound conflicting in a two body system, but in fact, it is a direct consequence of the flxed number of particles 
that characterizes the Canonical Ensemble approach. In Fig. [5] we plot the pressure in the homogeneous plateau 
as a function of both density parameters. In Fig [5] (a) we adopt the plateau density po. There, an important 
characteristic is apparent, for small cavities a/2 < R < a we obtain po = and Po — (see also Figs. [2] and [3]). 
The anomalous bi- valued behaviour of Po{po) is a consequence of the non monotonic behavior of the density in 
the homogeneous region as was shown in Fig.[2l In Fig.[5](b) we adopt the rough density p. For Poip) we flnd a 
simpler general dependence with a monotonic behaviour until p = 2/^(1) is reached. Note that the maximum 
attainable density is Pmax = 2/y(l/2). In both Figs. [5] (a) and (b) the pressure attains its maximum when the 
density plateau disappears at R = a and then pressure drops discontinuously to zero. Interestingly, both used 
density parameters {po,p} are usually utilized to describe the behaviour of macroscopic fluid systems. From 
an opposite point of view, we may concentrate in the external force and on contact properties. From the wall 
theorem |lip28j the total scalar force between the wall and the HS system in a pressure form is 

dR HQ) ■ A{Ry^ = PPw = Pc , (113) 

with Pc from Eq ((62l) . The general features of (3Pw{R) ior D — 2 and 3 can be seen in Fig. 1 of Ref. [l9]. There, 
the basic systematic behaviour with increasing dimensionality is apparent. As a consequence of the chain rule 
of the derivative and the sticky bond transformation the numerator of pc is the CI of one HS in a diminished 
volume. The allowed volume for the HS particle is the free volume available when the other HS is slicked to the 
surface. The asymptotic behaviour of Pw may be obtained from Eqs. I|53[ [58|) . For D — 2 and 3 it was studied 
in |i49]. In the caging limit, when R^ a/2, Pw diverges as 



^.£±lmv„-,r\ (115) 

where Vq — V{a/2) and the expressions in Eqs. (|114[ IllS)) are consistent to order minus one in {R— §) and 
(y/Vb — 1). The same power dependence of the compressibihty factor was found for the caging Hmit of N- 
HS systems under periodic boundary conditions [S]. A comparison between Eqs. I|107p and (|115p shows an 
interesting similarity at D = 3. 

Finally, we shall study the surface tension of the system. Due to the failure of the obtained pressure tensor 
we can not evaluate a microscopic expression for the surface tension, however we may adopt a macroscopic 
approach. Identifying the radius of the dividing interface with R we define 7 by 

l3n{R) = - PPo - P^P , (116) 

where at the left hand side appears the pressure difference AP — Pw — Pq. If we identify 7 with the surface 
tension, the Eq. (|116p is the original version of Laplace equation applied to a spherical substrate-fiuid interface 
[llj . Besides, for D = 1 or i? — *■ 00 we obtain the planar equilibrium condition pc — PPq. Figure [6] shows 7 
for different radii and several dimensions D > 1. Some general features are: it is a positive defined quantity; at 
large R the value of 7 goes to zero with an expected power law dependent on D; in the opposite, as i?/cr — + 1/2 
7 increases indefinitely as pc do. It can be observed a finite jump at R/a = 1 due to a discontinuity in Pq. 
Probably such discontinuities in Pq and 7 are unphysical artifacts of the adopted definition for Pati^)- For 
R > a we have 

^ 2[V ~ {1 ~ ciz))b - V + b] /Q , 

= 2&c(z)/Q, (117) 

c(z) = 1 - wd{z) - {2z)'^wd{1 - 22^) , (118) 
where < c(z) < 1 and the most relevant terms of c{z) are 



c{z)^{D-1)Bdz 



2(1? + l)-i 



£> - 3 
3(1) + 3) 



(119) 



where for Z) = 3 we obtain the somewhat surprising result c{z) = z3/4. To first non null order in density and 
curvature we obtain 



where first term does not show any curvature dependence. In Addition, to the same order of accuracy in density 
we may replace p ^ pq. We identify the first density-non-curvature term as (3^ fiat = a/2p^. To first non null 
order we obtain 

/cr\2 (15-3) (£1 + 1) , , 

7/7/..-l--y 24(I)V3) ' (1^1) 

where the right hand term becomes null for D — 3. Unfortunately 7 is not the surface tension but in fact it 
is an excess free work. A better proposal for the definition of the surface tension 7 is a refined version of the 
Laplace equation [2^ 

9fl7+j7 = Z\P, (122) 




R/a 

Figure 6 (color online) Function 7, its dependence on pore radius for several dimensions. The dashed line shows the 
discontinuity at R/a = 1. 



where AP = j^. Definition from Eq. I|122p separates an explicit curvature dependent term from AP jll]. The 
asymptotic behaviour of AP is essentially described in Eq. I|120p , from that we obtain for the surface tension 



/cr\2 (D- 1)(D + 1) , , 

where now, the right hand term becomes null at D = 1. We may study the same system from a different radii 
which is an interesting point of view principally for non sharp interfaces usually found in fluid droplets. Being 
i?" = i? + e and following Henderson [2^ we define, 



m /37(i?) = /37(i?, i?") + m"l{R. R") , (124) 



its solution is 



and then 



P^(R,R")^f3jiR)-^{D-l) ^—j }, (125) 

jiR,R")/j{R)-l = ^ ^ly + O^i^), (126) 

showing that first correction is second order in e/R and positive for D > 1. 

We expect that, several of the above expressions may be generalized for any number of particles TV with 
the addition of the overall A^-dependent correction factor [23] . Here we do not demonstrate but merely suggest 
that, the Eq. I|120p should be multiplied by 2(1 — N~^) at right hand side [26] to become valid to the same 
order in density and curvature. The same modification applies to 7/;at. However, this prefactor does not modify 
the ratios in Eqs. i|12H [T23l) . Therefore, for a large enough R we obtain for the A^-HS-HWSP system in D 
dimensions 

/?7/i,t~a(l-iV-i)p2, (^27) 

/cr\2 (D- l)(i:> + l) , , 

7/7/.*-l--y 24(I) + 3) ' ^'''^ 

both are new results even in the most interesting case of = 3. Furthermore, the Eq. I|127p represents the first 
term of the density power series of the surface tension for any hard wall cavity. In addition, Eq. (|128p applies 
too for the conjugate system of a hard wall spherical core surrounded by a HS fiuid. Finally, we are interested 
in establish some relations between the studied system and the open (grand canonical ensemble) system of HS 
in contact with a hard spherical wall. The low density limit of an inhomogeneous open system was studied 
by Bellemans and Sokolowski-Stecki [HHS] who found the surface tension virial series and apply it to the HS 
system in contact with a planar hard wall (cite I Ref. [8] and I, HI Ref. [45]) and contained in a HWSP (cite 
HI in Ref. [8]) in D = 3. In the HS inhomogeneous fiuid in a HWSP the 7 first power density coefficient and 
zero order in curvature is the same that appears when a planar wall is studied (from Area term in Eq. (|53|) ) 
/3g|" = Wi'^/AoiR) = ao = 6z3+i(27r)-i which may be found by taking ^ 00 at Eq. ifHTp . The first 
curvature correction of this should be (5^^ . Bellemans obtained that in three dimensions this constant is null, 
^nf^ — ^3^^ — Oi ^iid we obtain (5^' = for D > 2 (from the nonexistent term in Eq. l(53|) ). whereas the 

first non null curvature correction is S^j^^ (from R^~^ curvature term in Eq. l(53|) ). To our best knowledge, all 
this properties have been never studied or evaluated for D ^ 2 and 3 [8j[34 j l45 j l49j . and here we are showing the 
systematic dimensional dependence in terms of the second cluster integral coefficient in a higher dimensionality 
space (see Eq. ((55)) ). Concisely, Eq. (4.5) in III Ref. [8j for a D dimensional system must be written as 

r(l) 

ihfiat - 1 = + Oi{R-^) + 0{p) o{R-^) , 

an 

- \r) 233 (I? + 3) ■ ^ ' 

where now 7/;at is the surface tension of the open system in contact with a planar wall. An interesting fact is 
that first correction is then quadratic in curvature and zero order in density. Besides, for a HS fiuid in contact 
with a hard convex spherical wall (sometimes referred in the literature as a spherical cavity inside of the bulk 
fiuid) this first order curvature correction should be exactly the same. The term proportional to —R~^ in Eq. 
(|129p is 1/40 and 1/18 for Z) = 2 and 3 respectively and becomes greater than one for £> > 27. A consequence 
of Eq. I|129p is that the Tolman length of the athermal system of HS-HWSP scales as 5^°^ ~ p for low density. 
It seems that the study of the three particle system will provide the value of the proportionality constant. We 
may conclude this section by noting that a curvature correction proportional to a / R should exist at Eqs. (|12H 
I128P if we consider a somewhat more realistic potentials with soft repulsion and/or attractive well. 

Here, we have presented an study of the bulk properties of a system consisting in 2HS-HWSP. In such an 
analysis we followed deliberately a non-thermodynamic approach. Even when it may be unexpected, the direct 



derivation of properties such as the pressure and surface tension using a thermodynamic approach is not a 
simple task [50] . The full implementation of this path requires a careful evaluation of several problems related 
with the nonextensivity of the system. This subject will be analyzed in an incoming work. 



7 Final remarks 

The few body system consisting of two hard spheres in a spherical pore in any dimension has been studied in 
the framework of the canonical ensemble of statistical mechanics. It was showed that several properties of such 
an inhomogeneous spherical system can be exactly evaluated. Analytical exact expressions for the canonical 
partition function Qd-, density distribution function p(r), and pressure tensor Vir) were obtained. The iterative 
construction of these functions in terms of the same properties for dimensions lower than D was performed. We 
should emphasize that neither approximations nor power series truncation were done along these derivations. 
The study of the analytical properties of CI at the low density limit or large R value, and the highly confinement 
limit or final solid were analysed. We found that properties at both limits are correlated. This becomes clear 
for odd D values where the CI is a polynomial and low density limit involves high order R monomials, though 
caging limit properties relates with the degree of the zero of CI at i? = cr/2. We found that such zero is of order 

(i:> + 3)/2. 

Other systems which are closely related to the 2-HS-HWSP have been also analysed. The system of two 
sticky HS or rigid linear body into a spherical pore was tackled. The two HS into a spherical pore with a 
smaller internal and fixed hard core was studied with special emphasis in the limit of on surface confinement. 
Several exact relations between the three closely related systems were estabhshed by applying the sticky bond 
transformation which makes possible their unified study. It was examined the way in which properties of the 
three systems become strongly dependent on the low density regime and on the opposite caging limit. Several 
equations that relates the coefficients of these systems were obtained. 

The pressure tensor of the 2-HS-HWSP was investigated. The re-examination of the mechanical equilibrium 
condition for a system with a hard spherical boundary was done. New constraints between non ideal pressure 
tensor components and one body distribution function slope at contact with the curved wall was obtained. 
The analytical evaluation of one possible definition of the pressure were performed. The obtained expression 
disregards the former equihbrium condition and therefore the used pressure tensor definition must be considered 
incorrect or at most approximate. 

The EOS of 2-HS-HWSP system was studied. We have analytically evaluated the pressure-density relation 
and the surface tension. In connection with the open system of HS in a HWSP, our results are consistent 
with that obtained by Bellemans [8] the first correction in the surface tension due to the curvature in the 
confining surface is not of order R~^ in three dimension. We also show (see Eqs. (fSSj [55|) and l|129p ) that the 
first correction has order R^^ independently of the system's dimensionality and is zero order in density. In 
addition we determined the value of the coefficients corresponding to the first inhomogeneous and first non 
planar wall corrections for all dimensions, Eqs. (|54l [55]) . To the best of our knowledge, it is the first time that 
both coefficients are evaluated. 

The In-Out relation introduced in Ref. [1^ and described in Section [51 the stick transformation introduced 
in Section m and several results (e.g., Eqs. l(53l[66l [72]) and (|84l) ) suggest interesting Hnks with the mathematical 
theory of Convex Bodies also known as Integral Geometry [35]l37] that will be studied in future works. 
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A Appendix: Properties of WDiVi 

This appendix is devoted to the iterative construction of wu{y, The procedure to iteratively build wu{y, ^) 
is traced from Eqs. (f23t [31]) where the last one applies only for A < y < 1. We define the complementary 
function of ^^(y, Z\) (see Eq. l(3T|l ). w_D(t/, Z\) 

WD{y, + Af woix') - i (1 - Af wd{x") . (130) 



with x' = r' /2Ri = (y2 + A)/[y{l + A)] and x" ^ r" /2R2 = (y^ - A)/[y{l - A)]. From the definition of 
wniy, A) (pTI) . replacing wd{x) with the aid of Eq. l(23|) and rearranging terms we obtain 



WD{y,A) = (l + Z\2)t(;B_2(y,Zi) + 2ziw£,_2(2/,zi) 



i?2 



(131) 



Following a similar approach with Eq. I|130p and writing {-Ri, R2} in terms of y and A we obtain the iterative 
relations 



WDiy,A) = {1 + A^)wD-2{y, A) + 2AwD-2iy, A) - 
y^^^ {(l~y^){y' - A^)f''-'^'\BolD 



(132) 



WD{y, A) = 2A WD-2iy, ^) + (1 + A^)wD-2iy, ^) 



First functions of the series are 



^0(^,-4) 
wo{y,A) 



Ay-^ {{l~y^)iy' ^ A^)y''-'^^'2Bn/D 

w^,iy,A) = l/{l-A^) , 
w.,{y,A) = -A/il-A^) , 



1 

TT 

1 
TT 



arccos{ 
arccos{ 



jI±A_, 

y{i + A)' 

y^ + A ^ 



arccos{ 



arccos{ 



y{l + A)' 
wi{y,A) = 1 - y . 
wiiy,A) = -A{l-y) /y . 



jI_A_ 

y^-A 



(133) 

(134) 
(135) 

(136) 

(137) 

(138) 
(139) 



B Appendix: Properties of go (2) and 1*0(2) 

This appendix is devoted to deduct a few properties of functions qoiz) and ud{z). We begin summarizing four 
properties of Ix{a,b) (from [1] and [2] p. 944) 



Ix{a,b) = 1 - Ii-x{b,a) , 

140,6 da; = c ^x^I^ia.b) - c ^ / \ 4 a + c, 6 

y (a)7 (a + + c) 

4(a, 6) = Ix{a + 1, 6) + (a B{a, b))-^x''{l - xf , 

Ix{a,a) = 2"i/i_4(^_i/2)2(a, 1/2) , withx < 1/2 . 
The deduction of Eqs. (|4T| |44|) for g(z) follows from the definition on Eq. ifTS)) 



qoiz) = 2^ D / x-^^i WD(a;)dx , 
where < z < 1, which may be rearranged using Eqs. lfT9t I140p 

to(^) = 2^i? / t^-i/i_t2((i? + l)/2,l/2)dt, 

= 2^ (1 - z-°) - 2^-iL' / t^'^-^It{\/2,{D+\)l2)dt , 



(140) 
(141) 

(142) 
(143) 

(144) 



(145) 



integrating on through identity (|14ip . and using Eq. (|140p and definitions of Eqs. (fT9 j |4T1) 

qoiz) = 2^(1 - z^) - 2^ (i^/2/t(l/2, p + l)/2) - 2-^/t((i? + l)/2), (D + l)/2)) 

= 2^(1 - z^) - 2^ /i(l/2, (i^ + l)/2) + {2z f 1,2(1/2, {D + l)/2) + 

h{{D + l)/2, p + l)/2) - + l)/2, {D + l)/2) 

= -(2z)^/i_,2((i? + l)/2, 1/2) + {{D + l)/2, p + l)/2) 

= (2z)^u;c(2). (146) 

Interestingly, above expressions may transform to show a complete dependence on wdIz). Using the identity 
(|143p and the analytic extension of woiz) ifTOj) we find ud{z) — 2^^wd{z') with z' = 2z^ — 1, and 

qoiz) - 2-^wd{z') - {2z)^wd{z) , (147) 

therefore, qoiz) may be built in terms of the recurrence relation Eq. ((23|) . With the purpose of deduce the Eq. 
(|44|) we will find some identities not available in the literature of Refs. [2^1]. From (|140[ I142|) we obtain 



/i_,(a, b) = /i_,(a, 5 + 1) - (6B(a, b))-^x\l - a;)-^ , (148) 

/i_,(a, 6) = /i^,(a + 1,6) + (a S(a, 6))-ia;''(l - x)'^ , (149) 
applying both recurrence relations we find 

/i_,(a + 1,6+1)== /i_,(a, 6) + (a 6 B(a, 6)/(a + b))-^x\l - x)"(^- - x) , (150) 

a + 

/i_,(a + 1, a + 1) = /i_,(a, a) + ^^|^a::"(l " ^Ti^ - x) , (151) 

and then using Eq. l(4T|) 



z.,,(z) = .z,-2(.) + ^2((^ + i)/2) "''''(^ - -')'^-^^/'(^ - ■ (152) 

From Eq. (|4T1) the first functions ud{z) of the series may be obtained (see Eqs. (|42j [43]) ). Finally, we may 
mention that qoiz) could be also defined as the solution of the second order differential equation 

d^qo ~iD-l) z-'d.qo - 2^-'BnD {z\l - z^))^""''^^' , (153) 
with boundary conditions gn(0) = 1 and (/d(1) = 0. 
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